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Abstract: We study the viability of five-dimensional gauge theories as candidates for 
the origin of the Higgs field and its mechanism for spontaneous symmetry breaking. 
Within the framework of lattice field theory, we consider the simplest model of an SU( 2) 
gauge theory. We construct this theory on a five-dimensional orbifold which explicitly 
breaks the gauge symmetry to U( 1) at the fixed points of the orbifold. Using anisotropic 
gauge couplings, we find that this theory exhibits three distinct phases which we label 
as confined, Higgs and hybrid. Within the Higgs phase, close to the Higgs-hybrid phase 
transition, we find that the ratio of the Higgs to gauge boson masses takes Standard 
Model-like values. Precisely in this region of the phase diagram, we find dimensional 
reduction via localisation. 
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1 Introduction 

The Higgs sector of the Standard Model is an effective theory that seems to be con¬ 
sistent with observations, but at the same time hides a number of puzzles. One such 
puzzle is the origin of the potential responsible for Spontaneous Symmetry Breaking 
(SSB), which leads to the Brout-Englert-Higgs (BEH) mechanism [1-3]. Another is 
the mechanism behind the stability of the Higgs mass under quantum fluctuations; its 
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quadratic sensitivity to an ultra-violet cut-off is related to the so-called hierarchy prob¬ 
lem. A class of extensions to the Standard Model aimed at addressing these puzzles by 
the use of extra dimensions come under the heading of Gauge-Higgs Unification (GHU) 
[4-6]. In these models, the Higgs held originates from extra-dimensional components of 
the gauge held and gives rise to massive gauge bosons in the regular four dimensions. 
In this study, we consider the simplest case of one extra dimension where, due to the 
higher dimensional gauge invariance, the Higgs potential remains zero at tree level and 
is generated only through quantum effects [7]. Furthermore, the use of an orbifold 
geometry [8] can break the gauge group at the hxed points, allowing for a Standard 
Model-like Higgs held in the fundamental representation. 

In perturbation theory, where the extra dimension is compactihed and the hve- 
dimensional helds are expanded as a series of four-dimensional Kaluza-Klcin modes, 
the Higgs mass is independent of the cut-off (at one-loop) [9, 10]. However, we consider 
a pure gauge theory where it is known that the one-loop effective Higgs potential [11] 
does not develop a symmetry breaking minimum [12]. Although the Higgs mass typi¬ 
cally tends to be light [13], one way to circumvent this problem is to include fermionic 
degrees of freedom which induce SSB via the Hosotani mechanism [6, 14], Alterna¬ 
tively, by introducing an explicit cut-off in pure gauge theory, a one-loop potential with 
a symmetry breaking minimum can indeed be generated [15]; in the presence of this 
cut-off, the perturbative non-renormalisability of the five-dimensional theory can not be 
ignored, and the stability of the Higgs potential beyond one-loop has to be investigated. 

In the context of lattice held theory, five-dimensional GHU models were formulated 
in [16, 17] for an orbifold geometry. Subsequently, non-perturbative studies using Monte 
Carlo simulations of the path integral showed that SSB occurs in pure gauge theory 
[18, 19]. This observation was confirmed in [20, 21] via semi-analytic mean-held calcu¬ 
lations [22], An immediate concern here is the apparent disagreement between these 
results and Elitzur’s theorem [23], which states that a local gauge symmetry can not 
be spontaneously broken. However, their reconciliation was confirmed by the discovery 
of a global symmetry in [24]; it is the spontaneous breaking of this global symmetry 
which leads to the BEH mechanism in the context of non-perturbative GHU models 
[25]. Other previous explorations of this theory have focused on the case of a toroidal 
geometry which is known to posses confined and de-conhned phases separated by a 
hrst-order phase transition [26-30]; second-order transitions due to compactihcation 
were studied in [27, 31, 32], and the scalar spectrum was measured in [33]. 

In this paper we present the phase diagram and spectrum of the theory formulated 
using an orbifold geometry and an anisotropic lattice. We concentrate on the region 
where the the lattice spacing in the usual four dimensions, a 4 , is less than or equal to 
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the lattice spacing in the extra dimension, a^. Within this region, the theory exhibits 
three distinct phases, which we label as confined, Higgs and hybrid. We find a layered- 
like dimensional reduction in the hybrid phase, while in the Higgs phase, close to the 
Higgs-hybrid phase transition, we observe dimensional reduction on the boundaries via 
localisation. Significantly, within this region close to the Higgs-hybrid phase transition, 
we find that the spectrum resembles that of the Higgs sector. This theory is therefore 
a candidate to reproduce the Higgs sector of the Standard Model. 

2 Theoretical Framework 

This study will focus on a pure SU( 2) Yang-Mills theory in five dimensions. We label 
our continuum five-potential A M , where M — 0,1, 2, 3, 5, corresponding to the usual 
four-dimensions and an extra fifth dimension. We construct the theory using an orbifold 
geometry following [16, 19], where it was found that this theory exhibits SSB. 

2.1 Gauge Theory on the Lattice Orbifold 

In order to construct the theory on an orbifold, one first considers a gauge theory 
constructed on a five-dimensional periodic Euclidean lattice T x L 3 x 2 N 5l where T 
refers to the number of lattice points in the temporal direction, L to the number of 
lattice points in the spatial directions and 2Y 5 to the number of lattice points along 
the extra dimension. We denote the lattice spacing as a and label the coordinates of 
the lattice points via a set of integers n = {n^}. In what follows, // is used to index the 
standard four dimensions. The gauge field is defined as the set of parallel transporters 
(gauge links) { U(n,M ) e 517(2)}, identified by lattice coordinates and direction, that 
connect neighbouring lattice points. The orbifolding of the extra dimension is then 
achieved by a combination of a reflection 7 Z and group conjugation T g , hence leading 
us to the orbifolding condition 

(1 -KT g )U(n,M) = 0 , (2.1) 


which is essentially a Z 2 projection of the gauge links. The reflection operation acts on 
both the lattice points 



n = 7 Zn 

= (n M) -n 5 ) 

(2.2) 

and the gauge links 




7 ZU (n , / 1 ) 

= U (n, / 1 ), 

lZU(n, 5) = U\n — 5, 5) . 

(2.3) 

The group conjugation acts 

solely on the 

gauge links 



T g U(n, M) = 

g U (n, M) g- 1 , 

(2.4) 
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Figure 1. A sketch of the orbifold lattice and the gauge links: boundary 17(1) links are red, hybrid 
SU( 2) links sticking to the boundaries are magenta and bulk 517(2) links are blue. 


where g is a constant SU( 2 ) matrix such that g 2 is an element of the centre of SU( 2 ). 
Without loss of generality, we choose g = — ia 3 , and in accordance with gauge invariance 
the gauge group is broken down to U( 1) at the fixed points of the orbifold (ns = 0, N 5 ). 

The theory is now defined in the domain / = {n M , 0 < ns < N$} with volume 
T x L 3 x N 5 corresponding to Figure 1. We perform our study using a five-dimensional 
anisotropic Wilson action 


Qorb 



l 
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w tr {1 

fJt<v 


U^uin^)} + 7 ^tr{l - f/ M s(n M )} , 


(2.5) 


where w = 1/2 for plaquettes, [/ M „, living at the fixed points of the orbifold and 
w = 1 otherwise [ 16 ]. The anisotropy parameter allows for different lattice spacings 
a 4 and a 5 , where a 4 denotes the lattice spacing in the usual four dimensions and a 5 
denotes the lattice spacing in the extra dimension. In the classical limit 7 = a 4 /as and 
f3 = 4a 4 / g\, where 55 is the dimensionful continuum gauge coupling. I 11 what follows, 
we will use an equivalent pair of couplings /3 4 and ^5 which are related to the couplings 
in equation (2.5) via 

& = - ; @5 ~ fi'y ■ (2.6) 

7 

As depicted in Figure 1 , there are three types of gauge links: four-dimensional 17 ( 1 ) 
links residing at the fixed points of the orbifold, bulk SU( 2 ) links and extra-dimensional 
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links which have one end at a fixed point and the other in the bulk; these so-called 
hybrid links gauge transform as U —> at the left boundary (n 5 = 0) and 

U —> Q su ( 2 ' > Utt^ u< ' 1 ' ) at the right boundary (rz 5 = N 5 ). 

We note that the orbifold boundary conditions for the gauge field, when the ma¬ 
trix g is chosen to be the identity, are equivalent to open boundary conditions in the 
orbifolded direction, cf. [34], 

2.2 Spontaneous Symmetry Breaking on the Orbifold 

Global symmetries play a crucial role in a non-perturbative gauge invariant formulation 
of gauge theories. Elitzur’s theorem dictates that any physical effect associated with 
the breaking of a local symmetry must originate from the spontaneous breaking of a 
non-trivial global symmetry [23]. 

Within our set-up, the spontaneous breaking of the boundary U (1) gauge symmetry 
is governed by a so-called stick symmetry , S [24], This symmetry is defined by the 
combination S = Sl ■ Sr, where Sl is a symmetry defined on the left boundary via 

17(ra 5 = 0,5) —> g~ x U(n 5 = 0,5) (2.7) 

U (n 5 = 0, fj) —> g~ l U (n 5 = 0 , /i) g s . ( 2 . 8 ) 

The symmetry Sl, which is not a gauge or centre transformation, is generated by an 
element g s of the generalized Weyl group, Wsu( 2 )(U(l)) = Nsu( 2 ){U(l))/U(l), which is 
the coset of the normaliser of 17(1) in SU(2) divided by 17(1) [25]. It is a Z 2 symmetry 
(see Appendix A for a proof) obeying {g, g s } = 0 that acts only on hybrid and boundary 
links 1 . Sr can be defined on the right boundary in an equivalent fashion. 

By construction, the system has an additional symmetry, J u , defined as the reflec¬ 
tion about the centre of the extra dimension. It is important to note that S commutes 
with T ', and we only consider operators that respect T. Since it is odd under S, the 
order parameter for the spontaneous breaking of the 17(1) is 

Z? = Z k + T[Z k ], (2.9) 

where Z k is defined via equation (4.9). A non-trivial expectation value of Z[ signals 
SSB giving rise to a mass of the associated gauge boson. 

An important property of the orbifold, as shown in Appendix A, is that the trans¬ 
formation by the centre of 577(2), which is obtained by applying equations (2.7) and 
(2.8) twice, does not seem to carry any physical consequence. All operators associated 

1 When g is the identity there is no stick symmetry since it always commutes with g s . The stick 
symmetry appears only when the scalar left on the boundary by the orbifold projection is not in the 

adjoint representation. 
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with bosonic degrees of freedom that extend into the extra dimension are even under 
the centre symmetry, as they must contain an even number of hybrid links. Since there 
is no order parameter which is odd under this symmetry, finite temperature phase 
transitions (defined by the spontaneous breaking of the centre) cannot occur along the 
extra dimension of the orbifold. The direct consequence is that the masses associated 
with the orbifold operators can not be finite temperature (or Debye [35]) masses. 

Following the arguments presented in Appendix A, the stick symmetry appears 
to be valid only at finite lattice spacing. In the perturbative continuum limit 2 the 
symmetry, and hence SSB, disappears. This explains the absence of SSB in pure gauge 
perturbative models of GHU [12], As a consequence, the description of the Higgs 
mechanism in our model of non-perturbative GHU makes sense only as an effective 
theory. It also implies that if the bulk phase transition is of first-order then SSB 
persists in its vicinity. 

2.3 Monte Carlo Simulations 

Within the framework of lattice held theory, we approximate the path integral by means 
of Markov Chain Monte Carlo simulations. Our gauge held configurations are generated 
using a Hybrid-Overrelaxation (HOR) algorithm, which consists of one heatbath step 
followed by L/2 overrelaxation steps [36, 37], where L is the number of lattice points 
in the spatial direction. The heatbath algorithm is described in [38, 39] for the bulk 
577(2) links, and in [40, 41] for the boundary 1/(1) links. 

Simulating in the vicinity of a phase transition can cause critical slowing down 
of the updating algorithm, and can prolong the thermalisation process. In order to 
minimise this effect, the hrst part of our thermalisation process is augmented with 
extra overrelaxation steps, and in extreme cases, an additional “magnetic held-like” 
term similar to equation (4.7) is added to the action. It should be stressed that these 
techniques are employed only in the very hrst thermalisation updates. They are sub¬ 
sequently removed, and the system is allowed to “re-thermalise” before initiating any 
measurement of an observable. 

For our statistical error analysis, we use the package described in [42]; our quoted 
errors are the one-sigma statistical uncertainty coming from this package. 

3 Phase diagram and Dimensional Reduction 

An important step in understanding the properties of a theory is to determine the 
possible phases it can exhibit. We determine the location and order of the theory’s 

2 According to our results, this is the only continuum limit which exists, since the phase transition 
is first-order everywhere. 
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Figure 2. The left (right) panel shows example hystereses in the expectation value of the plaquette 
for fixed Pi = 1.94 on the four-dimensional boundary (bulk) layers. The red line corresponds to a hot 
start and the blue line to a cold start. 

phase transitions through the expectation value and susceptibility of plaquettes, U^m, 
determined after hot and cold starts; a hysteresis indicates a first-order phase transition, 
while a smooth cross-over must be accompanied by the appropriate volume scaling of 
the susceptibility in order to identify a second-order phase transition. We perform 
our study of the phase diagram on lattices of various volumes L J x 4 and find that, 
once L > 32, the location and order of the theory’s phase transitions remain constant. 
In Figure 2 we show example hystereses at fixed /3 4 = 1.94 on the four-dimensional 
boundary (n 5 = 0) and bulk (n 5 = 2) layers. In all layers we find two hystereses which 
correspond to the system’s two phase transitions, as shown in Figure 3. 

In this study, we concentrate on the region where the bare anisotropy 7 < 1, 
meaning that the ratio of lattice spacings < 24/05 < 1. Figure 3 shows the phase diagram 
for 04 G [1.5, 3.0] and j3 5 G [0, 2.0]. Within this region, we find no evidence of a second- 
order phase transition. Since we find only first-order phase transitions, the data points 
represent the location of the centre of the corresponding hysteresis. The bands plotted 
through the data show the width of the hystereses within that vicinity. 

As in the case of toroidal geometry [26] , we determine a first-order phase transition 
between the confined and de-confined phases represented by the blue line in the figure. 
We refer to this transition as bulk-driven since it is the SU{2) gauge held in the bulk 
that drives the system into a change of phase. We follow the norm and label a phase as 
confined (de-confined) when the expectation value of the Polyakov loop is zero (non¬ 
zero) in all directions. However, since massive Higgs and Z-like bosons have been shown 
to be a property of the de-confined phase [19], we refer to it as the Higgs phase , and 
delay the discussion of its properties until Section 3.1. 
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Figure 3. The phase diagram for N 5 = A in the region of the Higgs-hybrid phase transition. The 
points show the location of a first-order phase transition. The red and blue lines represent the width 
of the corresponding hystereses, while the dashed orange line represents 7 = 1 . 


We determine the location of a previously unobserved first-order phase transition, 
represented by the red line, which begins at /3 4 = 2.02, j3 5 — 0 and ends by merging 
into the bulk-driven phase transition at the triple point /3 4 = 1.9, /3 5 = 1.21. We 
refer to this transition as boundary-driven since it is the U( 1) boundary gauge fields 
that force the system into a change of phase. The four-dimensional U( 1) theory is 
known to exhibit a phase transition at f3u(i) = 1.01 [43], which corresponds directly 
to our value of /3 4 = 2.02 since our boundary plaquettes take w — 1/2 in the orbifold 
action, equation (2.5), thus forcing the location of the phase transition in terms of 
the four-dimensional coupling to double. This transition defines a third phase of our 
theory which we label as hybrid for reasons that will become clear when we discuss its 
properties in Section 3.2. 

One important observation is that the phase structure determined here agrees on 
a qualitative level with the one determined from mean-field calculations in [20, 44], 
Another important observation is that our phase diagram is reminiscent of the Abelian 
Higgs model in four dimensions for a Higgs field of charge q = 2 [45, 46]. This is ex¬ 
pected to be the theory to which the orbifold reduces in four dimensions [19]. In fact, 
the following sections will show that we find two regions at 7 < 1 where dimensional 
reduction occurs: the hybrid phase, and the region of the Higgs phase close to the 
Higgs-hybrid phase transition. This suggests a respective mapping of the hybrid and 












Dimension Yukawa 

Coulomb 

Confining 

4 c 0 Cl e “7 r 
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Cl 
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5 co Cl Ki< y r) 

co - 

£1 

r 2 
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Table 1 . Functional forms of four and five-dimensional Yukawa, Coulomb and confining potentials. 
mz corresponds to the mass of gauge boson associated with a Yukawa interaction, a corresponds to 
the string tension of a confining potential, Ci correspond to constants and K\ is a Bessel function of 
the second kind. 

Higgs phases of the orbifold to the Coulomb and Higgs phases of the Abelian Higgs 
model. Similarly, the confined phase of the orbifold reduces on the boundary, as —> 0, 
to the confined phase of the Abelian Higgs model. 

Extra-dimensional models are only phenomenologically viable if the extra dimen¬ 
sions are hidden at energy scales of the Standard Model. We are therefore only in¬ 
terested in regions of the phase space where dimensional reduction occurs to leave a 
Standard Model-like Higgs sector on the boundary. In order to determine the dimen¬ 
sionality of the system for a given set of parameters, we use the shape of the static 
potential ; we determine the potential, V(r), between static charges of the respective 
gauge group at a distance r located within layers orthogonal to the fifth dimension. 
The potential within a given layer is extracted from Wilson loops measured within that 
layer. In order to achieve an accurate determination of the potential, we replace the 
temporal links of the Wilson loop by their one-link integral [47] and employ the varia¬ 
tional procedure described in Section 4.2; in order to construct a large basis of Wilson 
loops we subject the spatial links to various levels of hyper-cubic smearing (HYP) [48]. 
We adapt this smearing to the orbifold and do not consider staples extending into 
the temporal or fifth dimension [49]. Once we have extracted V(r), we can determine 
its dimensionality and type by considering the global fitting forms given in Table 1. 
Subsequently, by considering r 2 F(r ) = r 2 V'(r), which is a renormalised quantity in 
four-dimensions, we can examine any deviation of the potential from the best fit in 
finer detail. 

3.1 The Higgs Phase 

We now turn our attention to the Higgs phase and study its physical properties in terms 
of the potential it exhibits. Throughout this work, due to their strikingly different 


9 









Figure 4. The left panel shows the potential, F(r), for the boundary (715 = 0) layer calculated on a 
volume 80 x 32 3 x 4 for = 2.1, ds = 1.075, which is located within the Higgs phase. The right panel 
shows the quantity r 2 F(r) computed for the same parameter set. The measurements are performed 
on 10, 000 gauge field configurations separated by 20 HOR steps. The fits described in the legend 
correspond to those in Table 1. 

qualitative behaviour, we will distinguish two distinct regions within the Higgs phase: 
7 = 1, and 7 < 1 close to the Higgs-hybrid phase transition. 

From measurements of the potential along the 7 = 1 line, we see no clear signs of 
dimensional reduction; the potential measured within the bulk layers strongly favours 
a 5-d Yukawa interpretation while on the boundary, the 4-d Yukawa fit is no more 
favoured than a 5-d Yukawa fit. However, as we move deeper into the 7 < 1 region, 
there is a tendency of the potential on the boundary layers to become more 4-d Yukawa- 
like, while in the bulk layers, the potential continues to strongly favour 5-d Yukawa fits. 
Significantly, once we reach low enough /3 5 , the potential on the boundary can only be 
described by a 4-d Yukawa fit. This is strong evidence that the method of dimensional 
reduction is that of localisation , that is, that the effective four-dimensional physics is 
localised to the boundaries. This mechanism is reminiscent of [50, 51]. 

In order to illustrate the statements above, the left panel of Figure 4 shows the 
potential, H(r), extracted from Wilson loops along the four-dimensional boundary 
(n 5 = 0) layer for /5 4 = 2.1, (3 5 = 1.075. The measurements are performed on 10,000 
gauge held configurations of volume 80 x 32 3 x 4 separated by 20 HOR steps. The 
potential clearly exhibits a 4-d Yukawa form. As shown in Table 2, the goodness of 
fit x 2 values per degree of freedom obtained for the other fits are at least an order of 
magnitude larger than that of the 4-d Yukawa. From the fit to the 4-d Yukawa form, 
we can extract the mass of the gauge boson that minimises the y 2 value. This turns 
out to be axrnz = 0.26. We also determine y' = (log(r 2 F(r)))' ~ a^mz [20] for the 
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Figure 5. The left panel shows the potential, V(r), for the bulk (715 = 2) layer calculated on a 
volume 80 x 32 3 x 4 for fa = 2.1, fa = 1.075, which is located within the Higgs phase. The right panel 
shows the quantity r 2 F(r) computed for the same parameter set. The measurements are performed 
on 10, 000 gauge field configurations separated by 20 HOR steps. The fits described in the legend 
correspond to those in Table 1. 

range of distances 3 < r < 7, and find a^nriz = 0.322(66). Both these determinations 
are consistent with the value extracted from the spectroscopic calculation described in 
Section 5, where we find a^rriz = 0.268(3). It is important to note that, the y 2 value of 
the 5-d Yukawa fit is minimised when mz = 0. In this case, the 5-d Yukawa fit trivially 
becomes of 5-d Coulomb type (see Table 1), and in order to distinguish the two in the 
figure, we constrain the mass of the 5-d Yukawa fit to be the one we extract from the 
spectroscopic calculation, namely a^mz = 0.268(3). 

The left panel of Figure 5 shows the potential extracted from Wilson loops along 
the bulk ( 77.5 = 2) layer for fa = 2.1, fa = 1.075. The measurements are performed on 
the same gauge field configurations described above for the boundary potential. In this 
case, the potential is clearly of 5-d Yukawa type; the y 2 values of the other fits are at 
least two orders of magnitude larger than that of the 5-d Yukawa. The Yukawa mass 
that minimises the y 2 value for the 5-d Yukawa fit is a^rriz = 0.25 which agrees with 
the mass we measure in our spectroscopic calculations; it seems that the SSB occurs 
in such a way that the mass of the gauge boson is a global property of the system. We 
remark that if mz is left as a free parameter in the 4-d Yukawa fit, it yields a value 
which is twice as large as the spectroscopic value; in Figure 5, we show the 4-d Yukawa 
fit with the latter mass. 

The right panels of Figures 4 and 5 show the quantity r 2 F(r ) which corresponds 
to the potentials shown in their respective left panels. In both cases, we take the 
best fit from the left panel and recast it into the form r 2 F(r). We then superimpose 
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Figure 6. The left panel shows the potential, V(r), on the boundary (715 = 0) layer calculated 
on a volume 80 x 32 3 x 4 for ,84 = 2.6, /I 5 = 0.6, which is located within the hybrid phase. The 
right panel shows the quantity r 2 F(r) computed for the same parameter set. The measurements are 
performed on 5,000 gauge field configurations separated by 5 HOR steps. The fits described in the 
legend correspond to those in Table 1. 


this fit onto the measured r 2 F{r). In both cases, we see only small deviations, if 
at all, of the measured r 2 F{r ) from its respective fit. Therefore, the correction to the 
boundary potential coming from the five-dimensional bulk is negligible, confirming that 
dimensional reduction has taken place on the boundary via localisation. 

3.2 The Hybrid Phase 

By sufficiently decreasing /3 5 and keeping /3 4 > 2.02, we enter a third phase of the theory. 
We label it as hybrid due to the fact that we observe four-dimensional de-confined 1/(1) 
theories on the boundaries, and four-dimensional confined SU(2) theories within the 
bulk layers; the layers in this phase appear to behave as independent four-dimensional 
slices that are very weakly coupled and hence, the physical content of a given layer is 
completely dominated by the gauge group of that layer. Furthermore, the Polyakov 
loop along the fifth dimension has zero expectation value, whereas the expectation 
value of Polyakov loops winding in the other four dimensions is zero in the bulk and 
non-zero on the boundary. This behaviour is strongly reminiscent of the layered phase 
described in [52], 

I 11 order to validate these claims, we show in the left panel of Figure 6 , the potential, 
P(r), extracted from Wilson loops measured along the boundary (n 5 = 0) layer for 
/?4 = 2.6, (3 5 = 0.6. From the figure it is clear that the potential can only be described 
by a 4-d Coulomb fit. As shown in Table 2 , the y 2 values for the other fits are two 
orders of magnitude larger than that of the 4-d Coulomb. On the other-hand, the left 
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Figure 7. The left panel shows the potential, V(r), within the bulk (n$ = 2) layer calculated on 
a volume 80 x 32 s x 4 for /?4 = 2.6, /?5 = 0.6, which is located within the hybrid phase. The right 
panel shows the quantity r 2 F{r) computed for the same parameter set, in units of the Sommer scale 
tq. The measurements are performed on 5,000 gauge field configurations separated by 5 HOR steps. 
The green lines (dashing) show the three-loop (two-loop) calculations of r 2 F{r ) in four-dimensional 
pure gauge 517(2) theory. The fits described in the legend correspond to those in Table 1. 


panel of Figure 7 shows the potential extracted from Wilson loops along the = 2 
bulk layer. In this case we find that only a confining string-like potential can describe 
the data; the other fits have y 2 values many orders of magnitude larger than that of 
the confining. Since we obtain confining potentials in the bulk layers, we can assume 
that the — C\jr term of the confining fit is the Liischer term, —(d — 2 ) 7 t/( 24 r ) [53, 54], 
giving us a probe of the dimensionality, d , of a given layer. For the potential shown in 
the left panel of Figure 7, we obtain a value of d — 3.97. Furthermore, we compute the 
slope, ci(r) = r 3 F'(r )/2 [55], which we find to plateaux at a value of —0.288(16). This 
is consistent with the expected value of — n/12 for a 4-d confining theory, emphasising 
the fact that the potential within a given layer of the bulk is indeed four-dimensional. 

When we attempt to compute the spectrum of the scalar channel within the hybrid 
phase, using the Higgs-type operators described in Section 4.1, we can not determine a 
particle mass. Instead, the signal coming from correlation functions almost immediately 
drop below machine precision signalling that the mass decouples from the theory at 
the scale of the hybrid phase; the bulk behaves like a two-colour version of QCD and 
we expect that the scale here is the Lambda parameter of 577(2). This expectation is 
justified since we can compute the Sommer scale r 0 [56], which is defined as a typical 
physical scale (0.5 fm) of QCD, from the potential shown in the left panel of Figure 7. 
We find that we reach the Sommer scale when the distance between two static charges 
is 8.87(5) in lattice units, meaning that the bulk layers within the hybrid phase are 
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Phase (Layer) 

4-d Yukawa 

5-d Yukawa 

4-d Coulomb 

5-d Coulomb 

Confining 

Higgs (Boundary) 

11.9 

112.7 

1.0 x 10 3 

138.0 

- 

Higgs (Bulk) 

918.7 

2.1 

4.0 x 10 3 

361.8 

- 

Hybrid (Boundary) 

208.9 

698.9 

9.9 

692.1 

- 

Hybrid (Bulk) 

1.6 x 10 5 

2.4 x 10 5 

9.7 x 10 4 

2.4 x 10 5 

2.2 


Table 2. X 2 values per degree of freedom for the fits to the static potentials shown in Figures 4, 5, 
6 and 7. 

sensitive to physics at the scale of two-colour QCD. Since we can not determine a scalar 
mass using our interpolating operators, and we know that glueball-like particles should 
be observable at the scale of two-colour QCD, we conclude that our operators do not 
have any significant overlap with glueball-like states and thus only probe the degrees 
of freedom relevant for Higgs-like particles, in agreement with [57]. We remark that 
the behaviour of the scalar and gauge boson masses, in the vicinity of the Higgs-hybrid 
phase transition, is consistent with the spectrum computed within the four-dimensional 
Abelian Higgs model [58]. 

By computing the quantity r 2 F(r), we can get a handle on the correction to the 
physics along the four-dimensional bulk layers coming from the extra dimension. The 
right panel of Figure 7 shows the quantity r 2 F(r ) at distances in units of the Sommer 
scale r 0 . One reaches perturbative scales as r/r 0 —* 0, and the green lines (dashing) 
correspond to three-loop (two-loop) calculations of r 2 F(r) in four-dimensional pure 
gauge SU (2) theory, see [27] . The fact that we see no deviation of our calculation from 
perturbative predictions down to very small values of r/r 0 tells us that the correction 
to the four-dimensional physics is negligible and occurs at energy scales much larger 
than the layers in the bulk are sensitive to. Furthermore, by recasting the confining fit 
shown in the left panel of Figure 7 into the form r 2 F(r) and superimposing it onto our 
calculated r 2 F(r), we see that there is no significant deviation from the 4-d confining 
fit. Therefore, it is clear that the correction to the bulk layers coming from the extra 
dimension is negligible. In a similar fashion, the right side of Figure 6 shows that the 
correction to the 4-d Coulomb boundary layers coming from the extra dimension is 
negligible. 

Significantly, on any four-dimensional layer within the hybrid phase, the 4-d Yukawa 
fit yields a zero mass for the gauge boson when the parameter mz is free. However, 
when mz = 0, the 4-Yukawa fit trivially becomes of 4-d Coulomb type. Therefore, in 
order to distinguish the two in Figure 6, we fix the mass in the 4-d Yukawa fit to be 
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a^iriz = 0.2, as it is a typical value we obtain for the Z boson mass in the lower region 
of the Higgs phase. The fact that mz = 0 minimises the y 2 of the 4-d Yukawa fit is 
strong evidence that there is no massive gauge boson and hence, no SSB in the hybrid 
phase. 


4 Spectroscopic Extraction 


We obtain spectral information from two-point correlation functions 


c l 3 (t) = - < o t (t))(o*m , ( 4 . 1 ) 


where 0*( 0) and Oi(t) are interpolating fields which correspond in the Hamiltonian 
formalism to operators that create and annihilate a state respectively. After inserting 
a complete set of eigenstates of the Hamiltonian, the two-point function becomes 


C«(i)=E 

n 


ryn* ryn 
Z C -E n t 

2E n 


(4.2) 


where E n is the energy of the n th state and the sum is over a discrete set of states due to 
the finite volume. The vacuum-state matrix elements Z = (n|(9j|0) can be interpreted 
as a measure of the contribution of an operator to a given state (see Section 4.2). 


4.1 Higgs and Z Boson Operators 

The main aim of this study is to reliably determine the low-lying spectrum of the 
J pc = 0 ++ (Higgs) and 1 (Z) channels within the Higgs phase. The quality of our 

spectral determination relies heavily on our ability to probe the relevant degrees of 
freedom of the theory. In order to achieve this in the Higgs channel, we follow [19], 
and begin by defining an extra-dimensional gauge potential. If the extra dimension is 
a discrete circle, it is given up to 0(a 3 ) by 

v M = ^(pM ~P f M) , ( 4 - 3 ) 

where p(n M ) is the Polyakov loop winding around the extra dimension. In order to con¬ 
struct this potential on the orbifold, S' 1 /Z 2 , the orbifolding condition (2.1) is imposed 
on the Polyakov loop, resulting in 

p(n^) = l(n / ,)gl f (n M )g f , (4.4) 

where l(n M ) is the Polyakov line extending across the extent of the extra dimension. 
Taking the potential, u(n M ), as our basic building block, we define the scalar field, 
/i(?V), as 

= H n „),g\ ■ (4.5) 
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o 

J PC 

5 

n 

0 ++ 

+ 

V 

0 ++ 

+ 

z 

1 “ 

- 

Z' 

1 “ 

- 


Table 3. The quantum numbers of our Higgs and Z boson interpolating operators. J PC denotes 
the total angular momentum, parity and charge conjugation in the continuum, while S refers to their 
parity under the stick symmetry. 


Since h{n^) is traceless, we define 

'H{n 0 ) = -jg J^tr [/i(n^)/i t (n M )] , (4.6) 

nk 

where the sum is over all spatial directions k, as our primary Higgs operator. From 
its construction, it is clear that it will act as a direct probe of the theory’s Higgs-like 
degrees of freedom. 

On the lattice, it is well-known that the use of the variational approach described 
in Section 4.2 leads to significantly improved spectral determinations. However, it 
requires the use of a reasonably large basis of interpolating operators, and hence, we 
add to our basis in the Higgs channel, a second type of operator 

tr [pM] • ( 4 - 7 ) 

We further increase our basis through the combined use of two smearing procedures. 
Firstly, we smear the gauge links i times using HYP smearing. We then construct 
Polyakov loops, pj(n M ), from the HYP smeared links U t , and apply j levels of an APE- 
like smearing to them 

Pii W = (1 - i) iP'/i) + C Q Ui ( n A n'JPUj-i)W U i ( n A*; rip) , ( 4 - 8 ) 

n 0 =n ' 0 
\n—n'\=a± 

where we take the constant c = 0.7, and the link n'J connects co-ordinate n /t to 

n'. Using these smeared Polyakov loops as the basic ingredients of equations (4.6) and 
(4.7), we arrive at our final operator basis for the Higgs channel, namely Vij and T-L^. 

Due to the nature of the theory’s SSB, we expect to determine a massive Z-like 
vector boson. Since it is the order parameter for SSB and possesses the correct quantum 
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numbers, the obvious choice of probe for this particle is the vector Polyakov loop. We 
define it as 

Z(n 0 ,k) = ^2 tr [gU (n M , + o 4 A:)17 t (n M , fc)a(n At )] , (4.9) 

rci,n2,n3 

where a(n M ) is the SU(2) projection of equation (4.5). As in the case of the Higgs 
channel, we can significantly improve our spectral determinations through the use of a 
large basis of interpolating operators. We therefore define a second type of operator, 
which we denote as Z' due to its importance in determining an excited Z'-like state, 
by 

Z'(n 0 ,k) =-^3 ^ tr \gU(n, k)l{n^ + aJi)U ] (n\ k)fi(nf)\ , (4.10) 

711,712,713 

where n = (n^n 5 = 0) and n' = (n /t , 77,5 = AT 5 ). Using the smeared Polyakov loops of 
equation (4.8), we arrive at our final operator basis for the Z channel, namely (Z k )ij 
and {Z' k )i. 

I 11 Table 3, we show the continuum quantum numbers of our interpolating opera¬ 
tors along with their parities under the stick symmetry S. Note that, in practice we 
construct our Higgs and Z boson operators at n 5 = 0, since we find that in general, 
(O) = (E[0}) holds. 

4.2 Analysis of Two-Point Correlation Functions 

In order to robustly determine spectral quantities, we employ a well-established vari¬ 
ational technique [59, 60]. Within a given symmetry channel, A, this amounts to 
constructing an N\ x N\ matrix of two-point correlation functions, and solving the 
generalised eigenvalue equation 

CijitX = A"(Mo)C^„K > (4- 11 ) 

where Cij(t ) is given by equation (4.1), and n = 1, 2,..., N\. Here, N\ refers to the 
number of interpolating operators used in a given channel. The solution procedure 
yields a set of eigenvalues {A n (t,t 0 )}, ordered such that A 1 (t, t 0 ) > A 2 (f,f 0 ) > > 

\ NA (t,t 0 ), and a corresponding set of eigenvectors { v for each t and reference time- 
slice to- At large enough values of t, the eigenvalues are proportional to and 

hence, the energy, E n , of state n can be extracted via 

04 m n = a^E n = In f -—-— ^ , (4-12) 

\X n (t + a iy t 0 )J 1 ; 

which is known as the effective mass of state n. 
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Figure 8. Effective masses of the ground states in the Higgs and Z boson channels along with the first 
excited state in the Z channel. The measurments are obtained using 5, 000 gauge field configurations, 
each separated by 5 HOR steps. The parameter set for this calculation is = 3.0, = 0.83, N 5 = 4, 

Nmggs = 12 and Nz = 8. 


In order to ensure that we determine reliable effective mass values, we take the 
following steps. Firstly, we construct an optimised operator basis in each symmetry 
channel. We achieve this by combining the operator overlap strategy described in 
[61], which uses vacuum-state matrix elements to distinguish between operators that 
significantly contribute to a given state from those that do not, with the absolute 
overlap strategy defined in [62], which gives an overall measure of how well an operator 
basis can reproduce a given eigenstate of the Hamiltonian. In general, we find that 
the addition of subsequent operators to our optimised bases do not significantly alter 
effective mass values of the low-lying spectrum. 

Secondly, we only determine effective mass values where a distinct plateaux between 
time-slices ti and tf is evident. By fitting the function c + e~ 5mt to the effective mass, 
we define ti to be the time-slice for which the exponential term of the fit is less than the 
one-sigma statistical uncertainty of the effective mass. We define tf to be the time-slice 
for which the value of the effective mass at tf is no longer within one-sigma of the fit 
function, or the one-sigma statistical uncertainty of the effective mass is twice that of 
ti . Furthermore, we only accept a plateaux that extends for at least four time-slices. 

In Figure 8, we show the results of applying the above procedure in both the Higgs 
and Z channels for the parameter set /3 4 = 3.0, (3$ = 0.83, N 5 = 4, and optimal bases 
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NHiggs = 12 and N z = 8 . The effective masses shown are obtained using 5, 000 gauge 
field configurations separated by 5 HOR steps, and are typical of all our calculations. It 
is clear that we reliably determine the ground state in both the Higgs and Z channels, 
and the first excited state in the Z channel. Although we see some evidence that an 
excited state in the Higgs channel may be present within the range 0.7 < a^E < 0.9, 
the effective mass does not meet the criteria described above and we do not show it in 
Figure 8 . 

5 Spectrum 

We now present the results of our spectroscopic calculations; preliminary results have 
already been presented in [44], Since we are interested in the properties of the theory 
where SSB occurs, we perform our spectroscopic calculations within the Higgs phase. 
Within the Standard Model, the ratio of the Higgs mass [63, 64] to that of the Z, 
p = mn/Triz ~ 1.38. In order for Gauge-Higgs Unification scenarios to be phenomeno¬ 
logically viable, they should achieve a similar value for a range of physically similar 
phase space parameters; an indication that this is achievable in our current SU{ 2 ) 
model would encourage further explorations into larger models that can account for all 
the degrees of freedom within the Higgs sector. 

In this section we fix N 5 = 4, and in Section 6 we take first steps in exploring the 
effect of varying the size of the extra dimension. Where instructive, we show our results 
in units of the energy scale 1 /R, where R is the radius of the extra dimension. However, 
it must be noted that, once an anisotropy is introduced into the couplings, the radius 
is only determined classically. In order to overcome this limitation, the renormalised 
anisotropy must be obtained, and work in this direction is already under way. 

Here we explore two distinct regions. The first being that of isotropic couplings, 
where /3 = f3 4 = ^ 5 , and the second being that of anisotropic couplings such that 7 < 1 . 
We show all of our spectroscopic results for these regions in Appendix C. 


7 = 1 

In Figure 9, we show our calculated Higgs and Z boson masses, in units of 1/R, for 
various values of the isotropic coupling, /5, and volumes, 80 x L 3 x 4; within all of 
our spectroscopic calculations, we fix the temporal extent T = 80. It is clear that we 
see no significant volume dependence for L > 32 and consequently, unless otherwise 
stated, we perform our calculations on volumes 80 x 32 3 x A 5 . Significantly, we find 
a^iTiz 7 ^ 0 as L —y 00 and hence, confirm the existence of SSB within the Higgs phase. 
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Figure 9. The Higgs (left panel) and Z (right panel) boson masses, in units of 1 /R, determined for 
various values of the isotropic coupling, /3, and volumes 80 x L 3 x 4. The black vertical line corresponds 
to the first-order phase transition at /3 ss 1.65 between the confined and Higgs phases. The orange 
horizontal line refers to the value of the one-loop Higgs mass from perturbation theory [9, 10]. 


The black vertical line represents the bulk-driven Higgs-confined phase transition. The 
orange horizontal line represents the Higgs mass calculated from a one-loop calculation 
in perturbation theory [9, 10], making it clear that the mechanism giving rise to the 
masses is highly non-perturbative and different from a finite temperature Debye mass, 
as discussed in Section 2.2. Clearly, for our range of /3 values, m H < rriz■ However, 
in the perturbative limit, where /3 — > oo, we expect the Higgs mass to approach its 
one-loop value and the Z to become massless, implying that a flip in the hierarchy is 
guaranteed. For our explored values of /3, we find p to be within the range 0.55 < p < 
0.75, suggesting that one would have to go to extremely large (and computationally 
impractical) values of /3 to achieve a value of p ~ 1.38, if at all possible 3 . We therefore 
rule out the possibility, at least for A r 5 = 4, of achieving a Standard Model-like hierarchy 
when 7 = 1. 


7 < 1 

We now turn our attention to the case of anisotropic couplings resulting in 7 < 1. 
The mean-field calculation presented in [21] suggest that our theory should produce a 

3 As we discussed at the end of Section 2.2, the continuum limit at (3 —)• oo is trivial. At finite 
values of /3 we only find first order phase transitions. Therefore the theory has to be treated as an 
effective theory at finite values of the cutoff. 
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Figure 10. The left (right) panel shows the ratio p = mn/mz for various values of (3$ and fixed 
(3 4 = 1.9 (2.1). The calculations are performed on various volumes, 80 x L 3 x 4, as described by 
the legend. The black vertical line on the left (right) panel shows the location of the Higgs-confined 
(Higgs-hybrid) phase transition, while the orange horizontal line shows the Standard Model value for 

P- 

Standard Model-like hierarchy close to the bulk-driven phase transition. Furthermore, 
it suggests that an excited Z' particle is present in the spectrum and that it is possible 
to construct renormalized trajectories along the phase diagram such that the physical 
quantities remain constant as the lattice spacing is varied. 

In order to explore this region, we start from fixed /3 values along 7 = 1 and 
determine the low-lying spectrum as we systematically decrease /3 5 . The left panel of 
Figure 10 shows the approach towards the bulk-driven phase transition for = 1.9, 
where it is important to note that, even close to the phase transition, we do not see 
any hint of finite volume effects. The black vertical line shows the location of the bulk- 
driven phase transition separating the Higgs from the confined phase, while the orange 
horizontal line shows the Standard Model value for p. It is clear that, in the vicinity 
of the phase transition, there is no qualitative change in the behaviour of p. 

On the other-hand, the right panel of the same figure shows the same procedure 
but for /?4 = 2.1. In this case, a major qualitative change in the behaviour of the 
theory is apparent, as significantly, p ~ 1.25 close to the phase transition. Note that 
in this case, in contrast to /3 4 = 1.9, the phase transition separates the Higgs from 
the hybrid phase. We find that on passing ^4 = 2.02, which is the location of the 
boundary-driven phase transition, the physical properties of the theory change such 
that the spectrum begins to resemble that of the Standard Model. We attribute this 
to the fact that beyond /7 4 = 2 . 02 , the (7(1) gauge field on the boundary naturally 
exhibits a de-confined behaviour, as discussed in Section 3. 
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Figure 11. The left (right) panel shows the ratio p = mn/mz for various values of /3s and fixed 
Pi = 2.7 (3.0). The calculations are performed on various volumes, 80 x L 3 x 4, as described by 
the legend. The black vertical line shows the location of the Higgs-hybrid phase transition, while the 
orange horizontal line shows the Standard Model value for p. 

The left panel of Figure 11 shows the approach towards the phase transition for 
fixed /?4 = 2.7. For /3 5 = 0.9, we obtain a value of p consistent with that of the Standard 
Model, p = 1.33(9), and in the right panel of the same figure we show that it is even 
possible for the theory to produce values of p significantly larger than the physical 
value of 1.38. Figure 12 shows a summary of our spectroscopic calculations close to the 
bulk-driven phase transition for 7 < 1 and N 5 = 4. As indicated by the legend, the 
different shadings of green correspond to three distinct regions in which we find p < 1 
(lightest shade), 1 < p < 1.3 (middle shade) and p > 1.3 (darkest shade). 

For various values of the couplings close to the bulk-driven phase transition, we 
robustly determine an excited Z' state. Within this region, we find that the Z' is 
roughly three times heavier than the Z. ffowever, we see evidence that, as we lower 7 , 
the separation between Z and Z' increases. This suggests that, in order to determine a 
more phenomenologically viable mass for Z', we must continue our exploration close to 
the Higgs-hybrid phase transition to larger values of /I 4 . We note that, although we do 
not robustly determine an excited Higgs state within our explored range of couplings, 
we see evidence that close to the Higgs-hybrid phase transition it has a mass around 
four times that of the Higgs. 
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Figure 12. Summary of our spectroscopic calculations for 7 < 1 and N 5 = 4 close to the Higgs- 
hybrid phase transition. The black points indicate the location of a spectroscopic calculation within the 
phase diagram. The green shading indicates the p = mu/mz value obtained for a given calculation. 
The lightest shade corresponds to p < 1, the middle shade to 1 < p < 1.3 and the darkest shade 
corresponds to p > 1.3. 

6 First Steps Towards Varying the Size of the Extra Dimen¬ 
sion 

As in all extra-dimensional models, it is crucial to determine how the physical properties 
of the theory depend on the size of the extra dimension. In this section, we take first 
steps towards answering this question by varying the parameter N 5 for a set of fixed 
couplings. Due to their qualitative differences, we again discuss the two regions 7 = 1 
and 7 < 1 separately. 


7 = 1 


The red triangles in Figure 13 shows the dependence of the ratio p on the parameter N 5 
for fixed isotropic coupling f3 = 3.0. As shown in Table 6, both the Higgs and Z boson 
masses remain relatively constant in lattice units, resulting in the constant behaviour 
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Figure 13. The ratio p for /3 = 3.0 (7 = 1, red triangles) and /?4 = 2.6, ^5 = 0.95 (7 ~ 0.6, blue 
circles) for various values of the extent of the extra dimension IV 5 . The orange line corresponds to the 
Standard Model value of p ss 1.38. The values were computed using 5, 000 gauge field configurations 
separated by 5 HOR steps. 

of p shown in the figure. We observe this behaviour for the range of isotropic couplings 
we explore and hence, it appears seemingly unlikely that the theory will produce a 
Standard Model-like p for any jV 5 while 7 = 1 . 


7 < 1 

We now turn our attention to the 7 < 1 regime, where we have observed that a 
Standard Model-like hierarchy is possible close to the Higgs-hybrid phase transition. 
We note that the location of this transition slightly moves into the Higgs phase when 
N 5 changes from four to six. Therefore, for the study of the N 5 dependence, we choose 
the couplings ^4 = 2.6 and (3$ = 0.95 as they are deep enough into the Higgs phase to 
avoid a change of phase as N 5 is increased. However, it should be noted that this shift 
seems to become negligible for higher N 5 values. The blue circles in Figure 13 show 
the dependence of the ratio p on the extent of the extra dimension N 5 . Table 7 shows 
that for increasing N 5 , the Higgs mass slowly decreases in lattice units. However, the 
mass of the Z boson rapidly decreases below the Higgs producing a Standard Model 
hierarchy for N$ = 8 ,10. This feature of the theory may become crucial in future 
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Figure 14. The masses, in units of 1 /R, of the Z (red triangles), Z' (magenta triangles) and Z" 
(purple circles) determined for a range of N 5 values at 84 = 2.6 and (3$ = 0.95. The magenta (purple) 
dashed line corresponds to the mass expected for the Z' ( Z") from a Kaluza-Klein scenario. 


work, where we attempt to construct lines of constant physics with a Standard Model 
hierarchy, as it opens up a much larger parameter space. 

An interesting observation is that the Z channel exhibits a Kaluza-Klcin-like be¬ 
haviour. Figure 14 shows the masses of the Z, Z' and Z" (where they could be robustly 
determined) for a range of N 5 values at 84 = 2.6 and /3 5 = 0.95. As shown in Table 7, 
the masses ~ I/A 5 in lattice units, producing the constant behaviour shown in the fig¬ 
ure. Furthermore, we observe that the masses of the Z' and Z" fit the excitation tower, 
mf K , expected from a Kaluza-Klein scenario. That is, that excitations mf^i appear at 
energy scales 1 /R higher than mf K . Although we see hints of an excited Higgs state 
which lies around 1/R above the ground state, we do not robustly determine its mass 
and hence do not show it in Figure 14. Work is under way in order to address this 
issue. 

7 Summary and Outlook 

We have explored the phase diagram and computed the spectrum of a five-dimensional 
pure SU(2) gauge theory formulated on a lattice orbifold using anisotropic gauge cou¬ 
plings. 
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In terms of the phase diagram, our results are summarised by Figure 3. We find 
a first-order bulk-driven phase transition separating confined and de-conhned phases. 
This line is analogous to the one present within a toroidal geometry. However, due to 
the existence of an extra global symmetry within the orbifold, a non-perturbative SSB 
mechanism is present, in accordance with Elitzur’s theorem. The de-conhned phase 
therefore becomes a Higgs phase due to the presence of massive gauge bosons. 

Since the orbifold breaks explicitly the gauge symmetry from SU( 2 ) to 17(1) at 
the boundaries, we find the appearance of a boundary-driven phase transition which 
separates the confined from a third phase of the theory. We label this phase as hybrid 
due to its properties; the boundary layers behave like a four-dimensional 17(1) theory 
in the Coulomb phase, while the four-dimensional bulk layers behave like a confining 
SU ( 2 ) theory. This phase exhibits a layered-like behaviour where each layer is (almost) 
independent of any other, and the physical content of a given layer is governed only by 
the gauge group of that layer. 

We have computed the spectrum of scalar (Higgs) and vector (Z) bosons in the 
isotropic 7 = 1 and anisotropic 7 < 1 regions. In the isotropic case, we find that 
the mass of the Z boson is always larger than that of the Higgs. However, in the 
anisotropic region close to the Higgs-hybrid phase transition, we find a Standard Model¬ 
like hierarchy of the Higgs and Z masses. This appears only to be possible once the 
boundary gauge held naturally de-confines, that is, once the four-dimensional coupling 
on the boundaries is larger than the critical value associated with the four-dimensional 
U( 1) phase transition. 

By computing the potential between two static charges in layers orthogonal to the 
extra dimension, we find that the system undergoes dimensional reduction precisely 
in the region close to the Higgs-hybrid phase transition where we find the Standard 
Model-like hierarchy; the four-dimensional Standard Model-like physics is localised on 
the boundaries, while the bulk remains Eve-dimensional. In the bulk, we End the po¬ 
tential to be of a 5-d Yukawa type, and the associated Yukawa mass to agree with the 
Z boson mass measured from our spectroscopic calculations. Significantly, the poten¬ 
tial on the boundaries exhibits a 4-d Yukawa form, where the associated Yukawa mass 
again matches that of the Z boson. 

The next steps in this project will address the issue of non-renormalisability of 
five-dimensional gauge theories. In the phase diagram we only see first-order phase 
transitions. On the other hand, we seem to need a finite cut-off for spontaneous sym¬ 
metry breaking to occur. We will therefore have to establish whether lines of constant 
physics (LCPs) exist for a finite range of lattice spacings. Furthermore, these LCPs 
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can be compared to those of the four-dimensional Abelian Higgs model, which is the 
model that this theory is expected to reduce to in four-dimensions. 

As a first step, we presented the dependence of the masses on the size of the extra 
dimension. The results appear favourable, and it seems that changing the size of the 
extra dimension will open up a much larger Standard Model-like parameter space for 
constructing LCPs. One requirement will be to determine the physical anisotropy, 
which is required to compute the physical size of the extra dimension. 

Work in the direction of the LCPs at 7 < 1, along with a study of the 7 > 1 regime, 
in which the theory is expected to compactify, is already under-way. 
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A Comments on the Stick Symmetry 

The stick symmetry Sl is defined by the transformations (2.7) and (2.8). Without 
loss of generality we can take g s = — icr 1 . Notice that while equation (2.8) is a Z 2 
transformation, equation ( 2 . 7 ) is a Z4 transformation. In particular the order 2 element 
of this Z4 takes a hybrid link into minus itself, yielding a transformation by a centre 
element. We denote the centre of SU(2) by Zf). The scalar and vector Polyakov loops 
defined in Section 4.1 are invariant under this Z!] (unlike on the periodic lattice). This 
specific part of the centre symmetry does not have any physical consequences, either of 
the finite temperature class or any other. Therefore the nontrivial part of the symmetry 
is obtained by taking the quotient Z4/Z!] which is isomorphic to Z 2 . Thus, the stick 
symmetry is indeed a Z 2 symmetry, identified as the generalized Weyl group [25]. A 
similar discussion applies to Sr. 
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We can attempt to follow the fate of the stick symmetry in the continuum limit. 
Assuming that the only continuum limit is the perturbative limit, we can write a 
hybrid link as U{n§ = 0,5) = e ia s( ulA \+ a2A l) where a§ is the lattice spacing along the 
fifth dimension. Then, defining </> = iA\ + A 2 and expanding in the lattice spacing and 
keeping only the leading order terms, we obtain from equation (2.7) 

1 —)• ia^cj) —> —1 

1 —>■ —ia 5 (j)* —ia 5 (J)* —>■ —1 (A.l) 

These transformations are consistent if A| = 0, that is if 0 = iA\. Had we chosen 
for g s the most general linear combination g s = —i (cos 6a 1 + sin(9a 2 ), we would have 
obtained the transformations corresponding to the linear combination (j) = {i cos 9 + 
sin0)A, where A is a real field. 

The transformation of the boundary links in equation (2.8) looks like a global [/(1) 
transformation so it has a well defined continuum limit. The two conclusions we can 
draw from the transformation of the hybrid links in equation (A.l) are that they are 
non-trivial only because the fluctuating degrees of freedom are group elements and that 
they interchange a field by a constant. The latter is consistent with the field taking a 
vacuum expectation value. The first conclusion on the other hand implies that in the 
naive, perturbative continuum limit one expects the stick symmetry together with SSB 
to disappear. Hence, the Higgs mechanism on the (pure gauge) orbifold lattice in the 
interior of the phase diagram seems to be tied to a non-zero lattice spacing which is in 
fact consistent with the presence of the first-order, bulk phase transition. The only way 
to circumvent this is if a second-order phase transition is found, where a non-trivial 
continuum limit could be taken. In such a case however one would have to possibly 
change the expansion of the link in terms of the continuum fields and the analysis would 
have to be repeated. Since at present we do not see any second-order phase transition, 
we will not discuss further this possibility. 

B Location of Phase Transitions 

In this appendix we give the location of the theory’s phase transitions (see Table 4) 
determined on L 4 x 4 lattices. The listed values refer to the centre of the hystere¬ 
sis, observed from the expectation value of the plaquette after hot and cold starts. 
The uncertainty associated with the values is given by half the observed width of the 
hysteresis. The plaquette expectation values were measured on 1000 gauge field config¬ 
urations, separated by 2 HOR steps. Note that volume checks have been performed and 
the phase transitions have been found to not show any dependence on L once L > 32. 
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Bulk-driven PT 


Bulk-driven PT 


L 

Pa 

(h 

32 

1.45 

1.95(10) 

32 

1.5 

1.85(10) 

32 

1.615(45) 

1.615(45) 

32 

1.70 

1.48(8) 

32 

1.80 

1.35(5) 

32 

1.90 

1.21(3) 

32 

1.92 

1 .20(1) 

32 

1.93 

1.19(1) 

32 

1.94 

1.188(13) 

32 

1.95 

1.17(1) 

32 

1.96 

1.1625(75) 

32 

1.97 

1.155(15) 

32 

1.98 

1.15(1) 

32 

1.99 

1.14(1) 

32 

2.0 

1.13(1) 

32 

2.01 

1.125(15) 

32 

2.02 

1 .11(1) 

32 

2.03 

1 .10(1) 

32 

2.04 

1 .10(1) 

32 

2.06 

1.08(1) 

32 

2.08 

1.065(5) 

40 

2.10 

1.05(1) 

40 

2.15 

1.015(5) 

40 

2.20 

0.985(5) 

40 

2.25 

0.955(5) 

40 

2.30 

0.93(1) 

40 

2.40 

0.890(5) 


L 

Pa 

P 5 

40 

2.50 

0.8625(75) 

40 

2.55 

0.8475(25) 

40 

2.6 

0.8387(38) 

40 

2.7 

0.8200(50) 

40 

2.8 

0.8050(50) 

40 

2.9 

0.7900(50) 

40 

3.0 

0.7800(50) 


Boundary-driven PT 


L 

Pa 

Ps 

32 

1.9 

1.21(3) 

32 

1.92 

1.175(15) 

32 

1.93 

1.155(15) 

32 

1.94 

1.137(13) 

32 

1.95 

1 .12(1) 

32 

1.96 

1.09(1) 

32 

1.97 

1.065(15) 

32 

1.98 

1.035(15) 

32 

1.99 

0.99(2) 

32 

2.0 

0.925(25) 

32 

2.01 

0.825(75) 

32 

2 .02(1) 

0.75 

32 

2 .02(1) 

0.5 

32 

2 .02(1) 

0.2 

32 

2 .02(1) 

0.01 


Table 4. The location of the first-order phase transitions for the IV 5 = 4 orbifold. Where indicated, 
the errors correspond to half of the size of the hysteresis. The triple-point is highlighted by the the 
light grey background. 
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C Tables of Masses 


In this appendix we tabulate, in lattice units, the masses obtained by the simulations 
described in Sections 5 and 6. Table 5 lists the masses measured on lattices where the 
extent of the extra dimension N 5 = 4, for different values of the couplings fa and fa. 
Tables 6 and 7 show the dependence on N 5 of masses measured at f3 = 3.0, 7 = 1, and 
fa = 2.6, fa = 0.95 respectively. 


fa 

fa 

Lattice volume 


CL^TYlz 

P 

a 4 m Z ' 

1.66 

1.66 

80x32 3 x4 

0.217(35) 

0.466(34) 

0.59(8) 

- 

1.66 

1.66 

80x48 3 x4 

0.246(16) 

0.418(18) 

0.59(8) 

- 

1.7 

1.7 

80x32 3 x4 

0.235(31) 

0.423(27) 

0.55(9) 

- 

1.75 

1.75 

80x32 3 x4 

0.241(31) 

0.343(56) 

0.70(12) 

- 

1.8 

1.8 

80x32 3 x4 

0.199(10) 

0.317(14) 

0.61(5) 

- 

1.9 

1.3 

64x48 3 x4 

0.204(25) 

0.433(5) 

0.47(6) 

- 

1.9 

1.3 

80x32 3 x4 

0.195(6) 

0.422(5) 

0.46(8) 

- 

1.9 

1.6 

80x32 3 x4 

0.239(26) 

0.361(34) 

0.66(15) 

- 

1.9 

1.9 

80x32 3 x4 

0.245(26) 

0.358(24) 

0 . 68 ( 8 ) 

- 

2.0 

2.0 

80x48 3 x4 

0.245(21) 

0.403(25) 

0.61(8) 

- 

2.1 

1.075 

80x32 3 x4 

0.327(8) 

0.268(3) 

1.22(4) 

- 

2.1 

1.1 

80x32 3 x4 

0.223(7) 

0.181(10) 

1.23(8) 

- 

2.1 

1.1 

80x48 3 x4 

0.215(19) 

0.183(21) 

1.19(12) 

- 

2.1 

1.2 

80x32 3 x4 

0.192(18) 

0.351(14) 

0.54(6) 

- 

2.1 

1.3 

80x32 3 x4 

0.243(7) 

0.405(23) 

0.60(5) 

- 

2.1 

1.3 

64x48 3 x4 

0.239(21) 

0.415(25) 

0.58(8) 

- 

2.1 

1.6 

80x32 3 x4 

0.237(22) 

0.337(23) 

0.71(8) 

0.550(12) 

2.2 

1.04 

80x32 3 x4 

0.231(14) 

0.242(19) 

0.96(8) 

- 

2.2 

1.09 

80x32 3 x4 

0 . 220 ( 11 ) 

0.299(24) 

0.75(7) 

0.49(5) 

2.3 

0.95 

80x32 3 x4 

0.215(5) 

0.187(2) 

1.15(3) 

- 

2.3 

0.97 

80x32 3 x4 

0.229(10) 

0.217(15) 

1.06(6) 

- 

2.3 

1.02 

80x32 3 x4 

0.239(11) 

0.236(12) 

1 . 01 ( 8 ) 

0.525(28) 

2.3 

1.3 

80x32 3 x4 

0.277(26) 

0.386(22) 

0.71(4) 

- 

2.4 

0.94 

80x32 3 x4 

0.214(17) 

0.204(8) 

1.05(9) 

- 

2.4 

0.99 

80x32 3 x4 

0.204(13) 

0.250(12) 

0.82(6) 

- 

2.5 

0.9 

80x32 3 x4 

0.197(8) 

0.184(5) 

1.07(6) 

- 

2.5 

0.95 

80x32 3 x4 

0.223(12) 

0.222(18) 

1.01(9) 

0.49(3) 

2.6 

0.85 

80x32 3 x4 

0.208(13) 

0.128(3) 

1.67(11) 

0.351(17) 
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At 

ft 

Lattice volume 

CL^TTIh 

CL^TflZ 

P 

CL^TYlz' 

2.6 

0.95 

80x32 3 x4 

0.238(9) 

0.240(8) 

0.99(7) 

0.58(5) 

2.7 

0.87 

80x32 3 x4 

0.233(23) 

0.175(8) 

1.34(12) 

0.546(25) 

2.7 

0.9 

80x32 3 x4 

0.267(15) 

0.201(8) 

1.33(9) 

0.37(5) 

2.7 

0.9 

80x48 3 x4 

0.250(10) 

0.208(16) 

1.20(9) 

0.550(17) 

2.7 

0.92 

80x32 3 x4 

0.26(3) 

0.22(1) 

1.15(13) 

- 

2.7 

0.97 

80x32 3 x4 

0.19(2) 

0.23(1) 

0.82(10) 

- 

2.7 

1.5 

80x32 3 x4 

0.221(7) 

0.274(15) 

0.80(5) 

0.424(5) 

2.7 

2.7 

80x32 3 x4 

0.212(18) 

0.316(12) 

0.67(7) 

- 

2.75 

0.86 

80x32 3 x4 

0.207(11) 

0.174(6) 

1.19(8) 

0.452(35) 

2.8 

0.91 

80x32 3 x4 

0.24(5) 

0.20(2) 

1.18(16) 

- 

2.8 

2.8 

80x32 3 x4 

0.236(12) 

0.356(18) 

0.66(5) 

- 

2.85 

0.862 

80x32 3 x4 

0.222(13) 

0.179(8) 

1.24(9) 

- 

2.85 

0.9 

80x32 3 x4 

0.216(11) 

0.207(8) 

1.04(6) 

- 

2.9 

0.85 

80x32 3 x4 

0.200(10) 

0.156(7) 

1.28(9) 

- 

2.9 

0.877 

80x32 3 x4 

0.26(2) 

0.216(4) 

1.19(1) 

- 

2.9 

0.9 

80x32 3 x4 

0.28(2) 

0.20(1) 

1.36(10) 

- 

2.9 

0.95 

80x32 3 x4 

0.274(12) 

0.256(8) 

1.07(6) 

- 

2.95 

0.892 

80x32 3 x4 

0.238(13) 

0.171(4) 

1.39(8) 

- 

2.95 

0.9 

80x32 3 x4 

0.22(3) 

0.21(2) 

1.06(18) 

- 

3.0 

0.83 

80x32 3 x4 

0.263(11) 

0.158(3) 

1.66(8) 

0.333(24) 

3.0 

0.87 

80x32 3 x4 

0.26(2) 

0.216(9) 

1.20(10) 

- 

3.0 

0.87 

80x48 3 x4 

0.249(14) 

0.193(7) 

1.29(8) 

0.491(25) 

3.0 

0.907 

80x32 3 x4 

0.269(15) 

0.241(6) 

1.12(6) 

- 

3.0 

0.92 

80x32 3 x4 

0.253(23) 

0.242(8) 

1.05(9) 

0.518(30) 

3.0 

0.97 

80x32 3 x4 

0.229(14) 

0.253(10) 

0.91(6) 

- 

3.0 

2.0 

80x32 3 x4 

0.203(14) 

0.323(19) 

0.63(6) 

0.583(11) 

3.0 

3.0 

80x32 3 x4 

0.251(7) 

0.335(11) 

0.74(3) 

- 

4.0 

4.0 

80x32 3 x4 

0.233(7) 

0.315(13) 

0.74(4) 

- 

5.0 

5.0 

80x48 3 x4 

0.192(16) 

0.297(15) 

0.65(6) 

- 


Table 5: Masses measured on JV 5 = 4 lattices. 
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N 5 

CI 4 TYIh 

a^rriz 

P 

4 

0.251(7) 

0.335(11) 

0.74(3) 

6 

0.232(12) 

0.289(13) 

0.80(5) 

8 

0.205(12) 

0.268(20) 

0.76(7) 

10 

0.219(24) 

0.297(11) 

0.73(9) 


Table 6. Masses determined at @4 = 3.0, = 3.0 for N 5 £ [4,10], on volumes 80 x 32 3 x N 5 . 


n 5 

a^rriH 

0,4771 z 

P 

fTlz 1 

O^TYlz" 

3 

0.284(12) 

0.325(4) 

0.88(4) 

- 

- 

4 

0.238(9) 

0.240(8) 

0.99(7) 

0.58(5) 

- 

6 

0.193(11) 

0.178(4) 

1.09(7) 

0.46(2) 

- 

8 

0.177(9) 

0.127(5) 

1.39(6) 

0.335(13) 

- 

10 

0.152(8) 

0.106(3) 

1.44(9) 

0.288(9) 

0.46(4) 


Table 7. Masses determined at 04 = 2.6, /? 5 = 0.95 for N 5 £ [3,10], on volumes 80 x 32 3 x N 5 . 
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